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Abstract 

We show that the quadratic interaction of fundamental and second harmon- 
ics in a bulk dispersive medium, combined with self-defocusing cubic nonlinear- 
ity, give rise to completely localized spatiotemporal solitons (vortex tori) with 
vorticity s — 1. There is no threshold necessary for the existence of these soli- 
tons. They are found to be stable if their energy exceeds a certain critical value, 
so that the stability domain occupies about 10 % of the existence region of the 
solitons. On the contrary to spatial vortex solitons in the same model, the spa- 
tiotemporal ones with s = 2 are never stable. These results might open the way 
for experimental observation of spinning three-dimensional solitons in optical 
media. 
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1 Introduction 



Solitons, i.e., self-trapped light beams or pulses that are supported by a balance 
between diffraction and/or dispersion and nonlinearity, are prominent objects in 
nonlinear optics Q. Optical spatiotemporal solitons (STS ) fl , alias superspikes 
[|| or light bullets Q| , were predicted in many works Q - |jl6| . They result from 
the simultaneous balance of diffraction and group- velocity dispersion (GVD) by 
self-focusing. Although they cannot be stable in the uniform self-focusing Kerr 
(x^) medium |B, stability can be achieved in saturable ||||, 10 1, quadratically 
nonlinear (x (2) ) [§, HI p), @ > 

and graded-index Kerr media []15| . STS can also 
be found in off-resonance two-level systems j^], in self-induced-transparency 
media flq| , as well as in engineered tandem structures incorporating quadrati- 
cally nonlinear slices |is|] . 

While a fully localized "light bullet" in three dimensions (3D) has not yet 
been found in an experiment, two-dimensional (2D) STS in a bulk \^ medium 
were observed in Rcf. p0|] . That work reported the formation of pulses in 
quadratic media, which overcome diffraction in one transverse spatial dimen- 
sion and GVD in the longitudinal direction. However, such experiments were 
performed by means of the tilted-pulse technique, which employs highly ellipti- 
cal beams; therefore, diffraction is negligible in the remaining transverse spatial 
dimension. 

Optical vortex solitons constitute another class of self-supporting objects, 
that have attracted much attention because of possible applications to the all- 
optical processing of information, or to guiding and trapping of atoms. The 
concepts of a multidimensional optical soliton and of an optical vortex may be 
combined, giving rise to spinning (vortex) solitons. Starting with the seminal 
works |2lj, both delocalized ("dark") and localized^ "bright" ) optical vortices 
were investigated in various 2D environments |^2|, |2^, Q |25|. In the 3D case, 
the bright spinning solitons take the shape of a torus ( "doughnut" ) [^6[ ^7) . 

For bright vortex solitons, stability is a major concern, as, unlike their zero- 
spin counterparts, the spinning solitons are apt to be destabilized by azimuthal 
perturbations. For 2D models with x^ 2 ' nonlinearities, an azimuthal instability 
was discovered by simulations |2g] and observed experimentally |2£| . As a result, 
a soliton with spin 1 splits into three or two fragments in the form of separating 
zero-spin soliton. Numerical simulations of the 3D spinning STS in the x^ 
model also demonstrates splitting into moving zero-spin solitons fj~7f . 

Nevertheless, the x^ 2 -* nonlinearity acting in combination with the self- defocusing 

(3) 

Kerr [x_ , where we use the subscript "minus" to stress the self- repulsion] non- 
linearity gives rise to stable spinning (ring-shaped) 2D solitons with spin s = 1 
and 2 j2j|. Models of this type for spatial [(2+l)-dimcnsional] solitons are well 
known |5(j| |3l|. The stability of the spinning solitons in the x^ '■ X-^ model 
may be realized as a result of competition between the self-focusing and self- 
defocusing nonlinearities. This understanding is further supported by the fact 
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that stable spinning solitons of the same type have also been found in another 
optical model displaying both focusing and defocusing nonlinearities, viz., the 
one based on the cubic-quintic (CQ) nonlinear Schrodinger (NLS) equation. In 
addition to optics, the same equations have been investigated in the contexts 
of Bose-Einstein condensates (BECs) j3^] and Langmuir waves in plasmas [ |33[ 
(however, in the former case, the quintic nonlinearity arises from three-body in- 
teractions, which also give rise to losses by recombination of BEC constituents 
into different species, thus making the quintic nonlinear coefficient a complex 
one). 

In the first direct simulations of 2D solitons with spin 1 in the CQ model, re- 
ported in the pioneer work |34], it was found that they are robust, provided that 
their energy is not too small |34fl. Later analysis, based on the computation of 
linear-stability eigenvalues, demonstrated that some of the spinning 2D solitons 
considered in Ref. Q are subject to a weak azimuthal instability. Nonetheless, 
in another part of their existence region, where they have a very large energy, 
the solitons with spin s = 1 and s = 2 were confirmed to be stable in the 2D 
CQ model |35| (see also Ref. for the stability investigation of the solitons 
with spin s = 1). Stable 2D vortex solitons in the CQ model can self-trap from 
Gaussian inputs with an embedded vorticity |37j. Notice that all the solitons 
with s > 3 have been demonstrated to be unstable in the CQ model p5j. 

A challenging issue is the search for physically relevant models in which stable 
3D spinning solitons exist. In fact, the only previously known model which could 
support stable 3D vortex solitons was the Skyrme model (see reviews |}8| ) . Very 
recently, we have found stable 3D spinning STS in the CQ model, which could 
again be construed as a result of the competition between self-focusing and self- 
defocusing j39j. Direct simulations of the 3D CQ model |f(| demonstrated that 
3D spinning solitons with moderate energies were unstable against azimuthal 
perturbations, while the ones with very large energies, i.e., broad "doughnuts" 
with a small hole in the center, were robust under propagation. However, a 
consistent stability analysis makes it necessary to compute eigenvalues of small 
perturbations. By calculating the instability growth rates, in Ref. [[[9) it was 
rigorously shown that sufficiently broad STS with spin s = 1 are stable, the 
stability region occupying ss 20% of their existence region, while all the STS 
with s > 2 are unstable. 

The aim of this paper is to show that the existence of stable spinning 3D 
solitons is a more generic fact, which is not limited to the CQ model considered 
in Ref. p9| . To this end, we will analyze the existence and stability of spinning 
STS solitons in the 3D version of the above-mentioned : X- model with the 
self-defocusing cubic term. In section 2, the model is formulated, and general 
results concerning the existence of 3D spinning STS in it, with different values 
of the spin, are displayed. Fundamental results for the stability of the spinning 
solitons, based on eigenvalues found from equations linearized around the soliton 
solutions, are presented in section 3. Direct simulations of the solitons' stability 
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within the framework of the full nonlinear equations are displayed in section 4, 
and section 5 concludes the work. 



2 The model and spinning solitons 

The scaled equations describing the reversible generation of the second harmonic 
(SH) from a single fundamental-frequency (FF) component u, in the presence of 
the self-defocusing cubic nonlinearity, dispersion and diffraction in the (3+1)- 
dimensional geometry, are well known |l2|, [l^, [2^, |3(], |3lj : 

du If d 2 u d 2 u d 2 u\ . .. ., . l2 , 

l az + 2 [up + w*- + m ) + u v - (H + 2H )u = °' 

dv If d 2 v d 2 v d 2 v\ „ , , , ., .,. 
l dZ + 4 [dx5 + dY~ 2 + a dT 2 ) ~ 0V + U 2(2|M| + |V| )V = °- (1) 

Here, T, X, Y and Z are the normalized reduced time, transverse spatial co- 
ordinates, and propagation distance, u and v are envelopes of the FF and SH 
fields, and f3 is a phase mismatch between the FF and SH waves. In particular, 
the variables used in Eqs. (0) are related to their counterparts (to be denoted 
by tildes) in Ref. [p3] as follows: u = 2u, v = w, V2(x, y) = (x,y). 

Equations ([!]) assume different GVD coefficients at the two harmonics, a 
being their ratio |Q, but neglects the Poynting- vector walkoff between the 
harmonics, and assumes that the temporal group- velocity mismatch between 
them |l4], [II], Q has been compensated. On the other hand, in the case a = 1 
the model possesses an additional spatiotemporal spherical symmetry |L2], 
Below, we will display results for the case a = 1, assuming that the group- 
velocity mismatch may be neglected in this case too. 

We look for stationary solutions to Eqs. (1) in the form u = U(r, T) exp {inZ + is0), 
v = V(r, T) exp [2 {inZ + is6)], where 9 is the polar angle in the plane (x,y), 
k is a wave number shift, and the integer s is the above-mentioned spin. The 
amplitudes U and V may be taken real, obeying the equations 

U pv , iw_^ v+ wy (2K+p)v+u ,_ 2(2U , +v2)v<ti) 
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dr 2 r dr r 2 dT 2 ) 



Dynamical equations (|l]) conserve the total energ 



Hamiltonian 

H = - 

2 



;y 

E = I I I (\u\ 2 + \v\ 2 )dXdYdT = E u + E v , (3) 



(\u x \ 2 + \u Y \ 2 + M 2 ) + ^(l^l 2 + \vy\ 2 + <j\vt\ 2 ) 
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+ [(3\v\ 2 ~(u* 2 v + u 2 v*) + (\u\ 4 + 4\u\ 2 \v\ 2 + \v\ i )] IdXdYdT, (4) 



momentum (equal to zero for the solutions considered), and longitudinal compo- 
nent of the orbital angular momentum L jl3| . The following relations between 
L, H and E for a stationary spinning STS follow from Eqs. L = sE, and 

H = -~kE + l -pE v -If J J (M 4 + 4M 2 M 2 + M 4 ) dXdYdT. (5) 

We have numerically found one-parameter families of stationary 3D spinning 
solitons which have the shape of a doughnut with a hole (supported by a phase 
dislocation) in the center. To this end, we solved numerically the coupled system 
of equations (||) using a standard band-matrix algorithm G3| to deal with the 
corresponding two-point boundary-value problem. We will display results for 
(7 = 1; however, we have also found that the STS exist for all a > 0, cf. Ref. 
[ fl3| , where nonspinning STS were studied in detail for a ^ 1. 

In Fig. 1 we summarize the output of extensive numerical calculations aimed 
to detect the domains of existence and stability of spinning STS. The continuous 
lines border the existence domain, and the dashed line constitutes a boundary 
between stable and unstable regions in the parameter plane (/?, k) . The way the 
stability boundary has been found will be explained in detail in the following 
sections. 

Shapes of three representative doughnut-like stable STS are plotted in Fig. 
2 for a fixed value of the net energy [see Eq. (S)], E — 12000. We see that, with 
the increase of the mismatch /3, the energy of the FF component of the spinning 
soliton increases, similar to the case of nonspinning solitons in pure media 

Below, we present systematic results which characterize spinning STS in the 
case of the zero phase matching, (3 — 0. In Fig. 3 we plot the curves n = k(E) 
and H = H(E) for both nonspinning and spinning STS in this case. The full and 
dashed lines in Fig. 3 correspond to stable and unstable branches according to 
results presented below. The s = solitons are stable according to the known 
Vakhitov-Kolokolov (VK) criterion, which states that the fundamental (s = 0) 
soliton branch undergoes a stability change at the point dEjdn = |48| . 

A feature shared by the nonspinning and spinning solitons, as it is evident in 
Fig. 3, is the absence of any finite threshold for their existence. This is a drastic 
difference from the recently studied STS in the CQ model, where well-defined 
thresholds were found for zero and nonzero values of the spin |20, p|] . 



3 Stability eigenvalues of the spinning solitons 

Complete understanding of the stability of solitons is provided by direct simu- 
lations of the evolution equations (see below) together with the analysis of Eqs. 
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(Q) linearized about the stationary spinning-soliton solution. In this section, we 
focus on the latter approach, seeking for perturbation eigenmodes in a general 
form, 

u{Z, r, T, 6) - U(r, T) exp [i(s6 + nZ)} 
= f(r, T) exp {\ n Z + i[(s + n)0 + kZ}} + g*(r, T) exp {X* n Z + i[(s - n)6 + kZ}} , (6) 

v(Z, r, T, 6) - V(r, T) exp [2i(s6 + kZ)} 
= p(r, T) exp {\ n Z + i[(2s + n)0 + 2kZ]} + q*{r, T) exp {\* n Z + i[(2s - n)6 + 2kZ]0?) 

where n > is an arbitrary integer azimuthal index of the perturbation, A n is 
the (complex) eigenvalue that needs to be found, and functions /, g and p, q 
obey equations 

Q T 2 



r (s + n) r f 

or 



-2 (U 2 + V 2 ) f - (U 2 - V) g - (2UV - U) p - 2UVq = 0,(8) 



-iKg + g 



i\ n p ■ 



-i\ n q + -. 



8 2 g 
dT 2 + 

-V 2 )g- 
d 2 p 
dr 2 
U 2 + V 2 ) 
d 2 q 

-4 (U 2 + V 2 ) 



n) 2 r 
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2 (U 2 
d 2 p 

-4 
d 2 q 



(2UV- 



0,(9) 



d 2 g 
— - 4 

dr 2 

(U 2 -V)f 

(2s + nfr 

or 

p - 2V 2 q - 2 {2UV -U)f- AUVg = 0(10) 



ng 

U)q- 2UVp 
- (2« + 0)p 



n) 2 r 



Or 



q - 2V 2 p - 2 (2UV - 



-{2n + 0)q 
U)g-AUVf 



0(H) 



Physical solutions must decay exponentially at r — ► oo. At r — > 0, f and g must 
vanish as rl s±n l, whereas p and q vanish as rl 2s±n l. 

To solve the above equations and find the eigenvalues, we used a known 
numerical procedure 126 
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which produces results presented in Fig. 4. The 
most persistent unstable eigenmode is found for value of the azimuthal index 
n = 2, for both s = 1 and s = 2. As is seen in Fig. 4, the instability of the 
soliton with s = 1, accounted for by ReA2, disappears with the increase of k 
at a stability-change point, K st ~ 0.04572, and the stability region extends up 



to K 



(3D) 
"-offset 



0.051, corresponding to the upper continuous line in Fig. 1, 
i.e., infinitely broad solitons (which implies that the vortex of the dark-soliton 
type [ pl| , that may be regarded as an infinitely broad spinning soliton, is stable 

too). The relative width of the stability region is ^k^^ — n st ^J /K^j^. ~ 0.1. 
However, no stability region has been found for the 3D solitons with s = 2, 
as well as in the 3D model of the CQ type, and in contrast to the 2D vortex 
solitons in both the '■ X— model with the competing quadratic and cubic 
nonlinearities (the same as considered here) p3[, and 2D CQ model ]3q|. 



G 



In the case when the spinning solitons are unstable, their instability is oscil- 
latory; the corresponding frequency, ImA [see Figs. 4 (c) and 4(d)] is found to be, 
generally, on the same order of magnitude as ReA at the maximum-instability 
point. In the stable region, k > K st , all the eigenvalues are purely imaginary. 
Oscillatory instabilities of solitons, characterized by complex eigenvalues of the 
corresponding non-self-adjoint linear operator, are typical to other conservative 
models of nonlinear optics |5(j, |5^] . 

4 Direct simulations 

The above results were checked against direct simulations of Eqs. (|lj), carried 
out by means of the Crank-Nicholson scheme. The corresponding system of 
nonlinear partial differential equations was solved by means of the Picard it- 
eration method |pq] , and the resulting linear system was handled by means of 
the Gauss-Seidel iterative scheme. For good convergence we needed, typically, 
five Picard iterations and fifteen Gauss-Seidel iterations. We employed a trans- 
verse grid having 121 x 121 x 91 points, and a typical longitudinal step size was 
AZ = 0.1. To avoid distortion of the instability development under the action 
of the periodicity imposed by the Cartesian computational mesh, we added ini- 
tial perturbations that were mimicking random fluctuations in a real system (cf. 
Ref. (10)). 

To illustrate the evolution of a stable 3D "bullet" generated by an input in 
the form of a completely localized Gaussian pulse with the energy Eq = 5986 
[see Eq. (||)], into which a vortex with s = 1 was embedded, in Fig. 5 we show 
the energies of its two components vs. Z. Robustness of the spinning STS is 
attested to by the fact that it can be generated from a Gaussian with a nested 
vortex, whose shape is far from the soliton's exact form. We see from Fig. 5 that 
there is a strong reshaping of the input Gaussian, which leads to a redistribution 
of the energy between the two components; some energy loss occurs, caused by 
emission of radiation in the course of the formation of the stable STS. Figure 
6 shows gray-scale contour plots of the intensity and phase distribution in the 
FF component, in both the input Gaussian with a nested vortex, and in the 
emerging spinning STS with the vorticity s = 1 at Z = 100, corresponding to 
the same case which was presented in Fig. 5. No further essential evolution of 
the soliton was observed in this case at Z > 100. 

Typical instabilities of the spinning STS with the spin s = 1 (in the case 
when it is unstable) and s — 2 are illustrated by Figs. 7 through 10. The 
azimuthal instability breaks the unstable spinning solitons into zero-spin ones, 
which fly out tangentially relative to the circular crest of the original soliton 
[similar to what is known about the instability-induced breakup of the (2+l)D 
spatial vortex solitons p8[]. Thus, the initial internal angular momentum (spin) 
of the doughnut-shaped spinning soliton is converted into the orbital momentum 
of the emerging nonspinning fragments. 
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Analyzing a large body of numerical results, we have concluded that the 
number of the emerging fragments is roughly equal to twice the original spin 
s. The dependence of the number of the fragments on the other parameters is 
fairly weak. 

It is noteworthy that, in all the cases displayed in Figs. 7 through 10 (and 
in many more cases not shown here), the number of the instability-generated 
fragments is exactly equal to the azimuthal index of the perturbation mode 
having the largest growth rate. Thus, the full nonlinear evolution of the unstable 
spinning solitons is in perfect agreement with the stability analysis based on the 
linearized equations, which was presented in the previous section. 

5 Conclusion 

In this paper, we have shown that stable bright spatiotemporal spinning solitons 
(vortex tori), which were recently found in the cubic-quintic model of a disper- 
sive optical medium with competing self-focusing and defocusing nonlinearities 
[ [39| , are also possible in a model based on the competition between the quadratic 
and self-defocusing cubic nonlinearities. The solitons are stable, provided that 
they are broad enough (so that the soliton's energy exceeds a certain critical 
value, or, in other words, the size of the internal hole is essentially smaller than 
the overall size of the soliton) . 

In fact, the model with the '■ X— (quadratic-cubic) nonlinearity may be 
realized easier in real optical media than the '■ X-^ (self-focusing-cubic - 
sclf-dcfocusing-quintic) one. Possibilities for the experimental implementation 
of the former model (chiefly, based on the quasi-phase-matching technique) were 
discussed in Refs. |2^, [34], [55| . Note that such optical media may be used equally 
well for the experimental generation of both the spatial (2+l)-dimensional soli- 
tons (vortex cylinders) considered in Ref. |^| and the 3D spatio-temporal 
spinning solitons (vortex tori) found in the present work. 

It is relevant to stress that the amplitude of a beam that can give rise to 
a stable spinning soliton should not be specifically large: as it is evident from 
Fig. 2, the necessary power is essentially the same as that which is necessary 
for the existence of a nonspinning soliton. The difference from the latter case 
is that the beam generating a stable spinning soliton must be broad (its cross 
section and temporal width should be large), i.e., its peculiarity is not a large 
power but rather large total energy. 

Similar to the cubic-quintic model, only spatiotemporal solitons with spin 
s = 1 may be stable in the present system, in contrast with the spatial (2+1)- 
dimensional solitons, which may be stable in the cases s = 1 and s = 2, in models 
of both types (on the other hand, a difference from the cubic-quintic case is that 
the existence of STS in the present model is not limited by any energy threshold) . 
These results suggest a conclusion that stable vortex solitons are generic objects, 
provided that the medium's nonlinearity contains competing elements and the 



8 



soliton's energy is large enough; in all the known models lacking the nonlinear 
competition, bright vortex solitons are subject to a strong azimuthal instability. 

Lastly, one can assume that, very generally speaking, the spinning soliton 
is not an absolutely stable object, but rather a metastable one. Indeed, the 
energy of the spinning soliton is larger than that of its zero-spin counterpart, 
hence it might be possible that a very strong initial perturbation will provoke its 
rearrangement into a zero-spin soliton, the angular moment being carried away 
with emitted radiation. In terms of this consideration, it appears that the s = 1 
and s = solitons are separated by extremely high potential barriers, which 
make the assumed process practically impossible. To illustrate this point, in 
Fig. 11 we show the cross sections of the s = 1 soliton which was very strongly 
perturbed at the initial point, Z = (the perturbation is w 30% of the soliton's 
amplitude), and the result of its evolution at the point Z = 200. For the same 
case, the comparison of the distributions of the intensity and phase inside the 
initial strongly perturbed soliton and the finally established one are shown in 
Fig. 12 (cf. Fig. 6). As is obvious from Figs. 11 and 12, the soliton was able 
to completely heal the damage, remaining a truly stable object. 
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Figure Captions 



Fig. 1. Domains of the existence and stability of spinning STS with spin s = 1. 
The upper continuous curve is the existence border, corresponding to infinitely 
broad (in fact, dark) solitons. 

Fig. 2. Typical shapes of stable STS with s = 1 for E = 12000: (a) 
= -0.1, (b) = 0, and (c) = 0.2. The labels FF and SH pertain to the 
fundamental-frequency and second-harmonic components of the soliton. 

Fig. 3. The propagation constant k (a) and Hamiltonian H (b) of the three- 
dimensional solitons, with different values of spin, vs. their energy E, in the 
case of zero phase mismatch, = 0. 

Fig. 4. The growth rate of perturbations, Re A, corresponding to different 
values of the azimuthal index n (indicated by labels near the curves) vs. the 
soliton's wave number k: (a) s = 1; (b) s = 2. The imaginary part of the 
stability eigenvalue, Im A, corresponding to different values of the azimuthal 
index n (indicated by labels near the curves) vs. the soliton's wave number k: 

(c) s — 1; (d) s — 2. Here and in the following plots, = 0. We stress that, in 
the case s = 1, the instability growth rate vanishes at the point the 
text, while in the case s = 2 the growth rate corresponding to n = 2 remains 
positive up to the border of the existence region of the solitons. This border is 
marked in all the panels by vertical arrows. 

Fig. 5. Evolution of the energy components E u and E v of the soliton with 
s = 1, as generated by an input configuration in the form of a Gaussian with a 
nested vortex. Here, the input total energy is E = 5986. 

Fig. 6. The formation of the soliton with spin s = 1 in the same case which 
corresponds to Fig. 5, shown in terms of the cross section of the fields at T = 0: 
(a) the intensity distribution in the initial Gaussian with a nested vortex; (b) 
its phase field; (c) the intensity distribution of the spinning soliton at Z = 100; 

(d) the phase field at Z = 100. 

Fig. 7. Isosurface plots illustrating the fragmentation of the s = 1 soliton 
with k = 0.01 into zero-spin ones as a result of the azimuthal instability: (a) 
Z = 0; (b) Z = 1000. 

Fig. 8. The same as in Fig. 7 in the case k = 0.032: (a) Z = 0; (b) 
Z = 1140. 

Fig. 9. The same as in Figs. 7 and 8 in the case of the s = 2 initial soliton 
with k = 0.015: (a) Z = 0; (b) Z = 900. 

Fig. 10. The same as in Fig. 9 in the case k = 0.04: (a) Z = 0; (b) 
Z = 2100. 

Fig. 11. Cross sections of an s = 1 soliton that was strongly perturbed 
at Z = 0, and the result of its evolution after having passed the propagation 
distance Z = 200. 

Fig. 12. The recovery of the soliton with spin s = 1 in the same case as 
in Fig. 11, shown in terms of the cross section of the fields at T = (cf. Fig. 
6): (a) the intensity distribution in the initial strongly perturbed soliton; (b) its 
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phase field; (c) the intensity distribution of the self-cleared soliton at Z 
(d) the phase field at Z = 100. 



15 



0.08 




-0.1 0.0 0.1 0.2 0.3 0.4 0.5 



mismatch 




energy energy 




wave number 



wave number 



5000 




1000 - 



20 40 60 80 100 
propagation coordinate 



This figure "Fig6.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/nlin/0206007vl 



This figure "Figl2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/nlin/0206007vl 



